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ABSTRACT 

Quenched chiral perturbation is used to study the decay constants for heavy-light mesons 
beyond the leading order in 1/M. The results are used to estimate the error in quenched lattice 
calculations of the decay constants. For the double ratio -Ri = (/d s //d)/(/s s //s) we find that 
the error is small — conservatively 5%, but likely smaller. We also find that quenching decreases 
the ratio /d„//d relative to the unquenched theory. 



PACS: 12.38. Gc, 12.39.Fe, 12.39.Hg 



1 Introduction 



The decay constants of the heavy-light mesons (B and D families) play a leading role in the physics of 
those mesons. In particular, B B mixing is proportional to f B . Thus, knowledge of fg is important 
for estimates of mixing rate and the extraction of the CKM angles. Consider the experimental 
quantity 

{AM/T) Bs 



To a good approximation, 
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However, in order for this to provide a meaningful determination of |V( S /Vtd|, the ratios fg B /f B and 
Bb s /Bb must be know precisely. They have been studied within heavy meson chiral perturbation 
theory (HChPT)[0, |, |, [| with the result that B B JB B is close to unity, while fgjfg is about 
1.2. The large deviation from unity of this latter quantity casts doubt on the reliability of HChPT. 
Because of this Grinstein|| has advocated using instead R\ x /d 8 //d, where 

Ri = ^/^- (3) 

SB 3D 

and JdJId is taken from experiment. 

The advantage of this approach is that the ratio R\ is sensitive only to corrections which violate 
both heavy quark and chiral symmetries. From the point of view of the heavy quark effective theory, 
this may be obvious, but it is nonetheless worth emphasising the idea behind it, which is simply 
that the long distance physics described by the chiral theory is insensitive to the actual mass of 
the heavy quark. Thus, one expects that while the chiral corrections may be badly behaved, their 
variation with respect to M (probed by Ri) should be under control. Indeed, this expectation was 
born-out in a recent calculation by Boyd and GrinsteinQ (BG) which found that i?i differs from 
unity by only —5%. 

The properties of heavy mesons have also been studied intently on the lattice. In particular, 
there have been many computations of fg and several of Bg • Thus the lattice provides an alternate 
way to determine R\. But of course it is necessary to know the error in these determinations. One of 
the systematic errors still present in most lattice calculations is that which arises from the use of the 
quenched (or valence) approximation, in which disconnected quark loops are neglected. Quenching 
alters both the short and long-distance properties^] of the theory. The short-distance the effects can, 
in general, be accounted for by appropriately adjusting the couplings of the theory, but the long- 
distance effects are more elusive. In a previous paper (|J), these effects were studied in heavy-light 
mesons by extending quenched ChPT||, ^[ |[ |j| to include heavy mesons at leading order in 1/M. 
It was found that quenching has a relatively minor effect on Bg and the Isgur-Wise function £ (w) , 
but a more pronounced effect on fg. However, the argument in the preceding paragraph suggests 
that the long-distance effects of quenching should largely independent of the heavy quark mass. 
Thus, one might expect that R\ and similar quantities will be relatively immune to the effects of 
quenching. 

In this paper we will attempt study this supposition in more detail. Our tool will again be 
quenched ChPT. Building on the work of BG, we will extend this to include 1/M corrections and 
use it to examine the 1/M dependence of the quenched corrections. By comparing these corrections 
to those found in unquenched ChPT by BG, one obtains an estimate of the error due to the use of 
the quenched approximation. Of course, the conclusions drawn from this approach are only reliable 
to the extent that ChPT accurately describes the unquenched physics. The remainder of this paper 

1 Here long and short distances are defined relative to the QCD scale. 
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is organized as follows. In section 2, we briefly review quenched ChPT for heavy mesons. In section 
3, we compute loop corrections to the heavy meson decay constants. In section 4 these results are 
investigated numerically. In section 5 we conclude. An appendix collects some results on the one- 
loop counter-terms. For the sake of brevity we will follow the notation of BG as much as possible 
and refer the reader to their paper for omitted details. 



2 Quenched Chiral Perturbation Theory and the Inclusion 
of Heavy Mesons 

Quenched chiral perturbation theory has been developed in refs. ^ and the extension to heavy 
mesons fields has recently been discussed in Ref. |5|. Consequently we will restrict the presentation 
here to only that which is necessary to fix notation. 

On way to implement quenching |l0| is to introduce bosonic "ghost" quarks to cancel the func- 
tional determinant which arises from the integral over the fermion fields. The inclusion of these 
extra particles enlarges the symmetry of the theory. Thus, quenched ChPT is obtained from or- 
dinary unqucnched ChPT by enlarging the symmetry group SU(3)l X SU(3)r to the semi-direct 
product (SU(3\3)l X SU(3\3)b)@U(1). Elements of the graded symmetry group are represented by 
supermatrices (in block form) 



C D 

where A and D are matrices composed of even (commuting) elements and B and C are composed 
of odd (anti-commuting) elements. Group invariants are formed using the super trace str and super 
determinant sdet, defined as 

str(f/) = tr(A) - tr(£>), (5) 

sdet(J7) = exp(str log (£/)) = det {A- BD^C)/ det(D). (6) 

To accommodate the larger symmetry, the meson matrix is extended to a supermatrix: 



,X 

where ~ QQ, X ~ QQ an d 4> ~ QQ and cf> is ordinary meson matrix 

J r^° + J r 



+ ^ K ° I. (8) 



k- k -y/iv 

Note that \ and x are fermionic fields, while cf> and (j> are bosonic. 
The Lagrangian of quenched ChPT is given by 

£2 2 

£Qx = Y [Md^d^) + 4 M0 str(M + )] + |^<M^o - (9) 
Here £ = £ 2 , £ = e^^ 1 (the normalization is such that f n = 128 MeV) while 

$ = -L str cj^-L^'-rf), (10) 

/ ra u \ 
M = m d , (12) 
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and M± = i(£ t .M£ 1 '±£.M£). The chief difference between the quenched and unquenched theories is 
the presence of the terms involving $ - In the unquenched theory they can be neglected because they 
describe the dynamics of the rj' meson, which decouples from the theory. But quenching prevents 
the rf from becoming heavy and decoupling, so these terms must be retained in quenched ChPT. 

The propagators that are derived from this Lagrangian are the ordinary ones, except for the 
flavor-neutral mesons, where the non-decoupling of $o leads to a curious double-pole structure. It 
is convenient to adopt a basis U a for the those mesons corresponding to mm, dd and so on, including 
the ghost quark counterparts. Then the propagator takes the form 

SjjCj (-aoP 2 + mg)/3 

ij p 2 - Mf + (p2 - MfW ~ Mf) { ' 

where e = (1, 1, 1, —1, —1, —1) and Mf — 2/iomi. It is conventional to treat the second term in the 
propagator as a new vertex, the so-called hairpin, with the rule that it can be inserted only once 
on a given line. Large N c arguments suggest that the kinetic coupling ao is small; this is supported 
by an analysis of 77 — 77' mixing. Consequently, we follow the usual practice and set it to zero in the 
sequel. 

In heavy quark effective theory Jll], [l2| [l3| the B and B* fields (for convenience we refer to all 
heavy mesons as B's) are grouped into the 4x4 matrix H a which conveniently encodes the heavy 
quark spin symmetry: 

H a = \(l + mTl,-B al5 ], (14) 
H a = 7 °^7°- (15) 

Here is the four-velocity of the heavy meson, the index a runs over the light quark flavors, m, 
d, s and the subscript "D" indicates that the trace is taken only over Dirac indices. Henceforth 
we will drop explicit reference to the heavy meson velocity. The quenched heavy mesons can be 
incorporated into this framework by adding to H extra fields B and B* derived from the heavy 
fields B and B* by replacing the light quark with a ghost quark, so that a also runs over the ghost 
flavors. 

BG (see also Ref. |l4|) have given the chiral heavy quark Lagrangian to order 1 JM (throughout 
the paper M will refer to the spin-averaged meson mass, M = \{Mb + 3Mb»)). To formulate the 
quenched version one must also include in the Lagrangian vertices which couple $0 to H . Symmetry 
requires that this coupling occur through str(A M ), which no longer vanishes. Adding these vertices 
and some extra counter-terms, but including only those terms which actually contribute at 0(1/M), 
one finds 

C = -(1 + jj) tv D [H a iv ■ D ba H b ] + g tr D \H a a^v ■ D ba H b a^} 
+ (g + ||) trc \H a H b Abalz] + H tr D [H a AbaJ 5 H b ] 
+ (7 + jj) tru [H a H a%ll5 ] str(^) + g tr D [H a7 ^ 5 H a ] str^) 
+ Ai tr D \H a H b ] (M + ) ba + ^| tr D [H a a^H a a^] + ^ tr D [H a a^H b <j^] {M+) ba 

jyi iVi 

+ (fcio + -77) tru [Haiv ■ D bc H b ] {M + ) ca + ^ tr D [~H a a^v ■ D bc H b a„ v ] {M + ) ca . (16) 
M M 

Light mesons enter this Lagrangian through the quantities: 

V M - l&^+^d^), (17) 
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A, = \{^-^d^)=-- f d^ + 0{^), (18) 

The couplings shift the masses of the heavy mesons: the B*-B mass splitting A is given by 
A = — 8A2/M, while Ai produces an SU(3) violating mass shift 8 q = 2\\m q and A3 violates both 
symmetries (though to the order we are working it only contributes at tree level). Taking these 
shifts into account the B q and B* propagators become 2 ^ v . k+ l A _ s ^ and ^frzi~^n > res P ec tively. 
The propagators of the ghost mesons are the same as their real counterparts. It is convenient to 
follow BG in adopting a shorthand notation for the Lagrangian (with the understanding that the 
propagators are the shifted one): 

£ = - tr D \H a iv ■ D ba H b ] + ~ giiH tr D [~H a H b 4 ba -f 5 ] + 7 Fff tr D [~H a H a 7^75] str(A"), (19) 

where 

[ 9b* =9+ 17 (.91 + 92) for B*B* coupling, 
9 = \ | (20) 

{ 9b =9+^(91-92) for B*B coupling, 

and similarly for 7. In the sequel it will be useful to introduce two more couplings, 

, (51+92/3) „ , (71+72/3) , . 

9b* =9 + j- , 1b*=1+ -j ■ (21) 

All of these couplings should be regarded as shorthand. That is, when performing arithmetic with 
these couplings, only terms to 0(1/M) are to be retained. Finally, we record the expression for the 
weak current] 2 ] (again with counter-terms) 

J£ = + £0 iXD^Ht^J + tr D [r^ la H b ^ ba } + ^ tr D [\r*, ^iD^H^l 

+ *a(^io + ^)tT D [r^H c C b }(M + ) ba + ia^tr D [^ la H c C b }(M + ) ba , (22) 
where r M = = 7^(1 - 7 5 )/2. 



3 Loop Corrections 

In ChPT, loop corrections generate terms which are non-analytic in the mass parameters of the 
theory. These terms are uniquely predictions because they cannot be canceled by higher-order 
counter-terms. The loops also generate divergences and finite terms which are analytic; these are 
absorbed into the counter-terms. In principle, they can be determined by fitting to experimental 
data, but due to the paucity of data in the heavy-meson sector, we will at times be forced to simply 
neglect them and they will not be shown in our formulas. We will later discuss the uncertainty 
arising from this approximation, but fortunately the error in R\, which is our main interest, should 
be free of counter-terms. 

As noted in Ref. ||, the loop structure of the quenched theory is rather odd: the theory splits 
into three copies of a one-flavor theory, distinguished only by the light quark mass. Because of this, 
all results for Bd or B u (henceforth we will refer to these simply as B) mesons apply for B s mesons 
after an appropriate relabeling. 

The loop integrals encountered arc cither identical to the 1\ and J of BG, or may be obtained 
from them by differentiation with respect to m 2 , which will be denoted with a prime. More explicitly, 

2 Here and in the Lagrangian |l6| several terms proportional to str(A1 + ) have been since they do not contribute, 
but their presence can be inferred from the gaps in the numbering of the coefficients. 
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Figure 1: The diagrams which contribute to the heavy meson self energy. Solid lines represent 
heavy mesons (J5, B*), dashed lines represent light mesons and the cross represents an insertion of 
the "hairpin" vertex. 

we have J(m, A) = AJi(m, A) and 



Ji(m) 



m 2 ln(m 2 //U 2 ) , 



Ji(m,A) = (-m 2 + |a 2 ) ln(m 2 /M 2 ) + ~ m 2 )F(m/A) , 



(23) 
(24) 



F(x) = 



a/1 — x 2 tanh 1 \/l — x 2 , x < 1, 
— \Jx 2 — 1 tan" 1 \J x 2 — 1 , x > 1. 



(25) 



The graphs which contribute to the self-energy are shown in Fig. [T]. In the diagram of Fig. [TJa, the 
ghost mesons will cancel the contribution from the real mesons unless one of the vertices involves the 
singlet field. Combining this with the contribution of the hairpin vertex Fig. [j}}, we obtain (recall 
that we are taking ao to vanish) 



i Tib(v • k) 



6i 



16tt 2 / 2 



A 



A 



2g B Zf B J{M d , - + S d -vk)+ g 2 »J>(M d , - 



v ■ k) 



(26) 



The terms not shown are analytic in M d and are absorbed in the definition of the counter-terms. 
For the B*, we similarly obtain 



«5«--*> = iff? 



u • £;) 



+ 2g B ^ B J(M d , + <J d - v • fc) + 4g B ^ B *J(M d , ^+S d -vk) + 



,(27) 



where we have neglected a term proportional to v^v 1 ' . Note that the 5 d dependence will vanish when 
the on-shell conditions are invoked. 

3.1 Wavefunction Renormalization and Decay Constants 

The wavefunction renormalization constants are obtained by differentiating the self-energy with 
respect to 2v ■ k and evaluating on-shell. BG made no assumptions about the magnitude of ^ L 
and thus retained the function F. In lattice simulations, where both and A are adjustable 
parameters, similar restraint is warranted and so we will also retain F. In terms of the shorthand 
couplings introduced earlier, the renormalization constants can be written compactly as 



Zr = 1 



o~2 ™o/3 
6 9b 



2 In A 2 T7f M d 



16tt 2 / 2 



,2 



log— f 



6g 2 



1/3 A 2 F(^) 
16tt 2 / 2 Mj - A 2 



24 37 



16tt 2 / 2 
A 2 



16tt 2 / 2 ^ A h 



(28) 
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(a) (b) 

Figure 2: The tadpole correction to the weak current vertex. The box represents an insertion of 
the weak current. 



j B « = 1 + -sg B 



ml/3 . Mj A A Mj A/2 

log -f + 6 5B . 7B* -^72 log 



., „z 2 /3 A 2 F(^) A 2 M d 



The above expressions have been simplified by discarding terms that are formally of higher order 
in the chiral and large mass expansions. However, it is important to retain those terms which are 
needed to make the chiral limit consistent. 

Loop corrections to the left-handed current vertex arise from the diagrams of Fig. 0. As was 
shown by BG, the diagrams of Fig. [2|a do not contribute at 0(1/ M). The remaining tadpole graph 
Fig. |b yields for the B 

iav»„ , pi + 2p 2 mg/3 Mj 
M J 16tt 2 / 2 g tj? 



and for the B* 



(l + ^^)^log^- (3D 



, Pr -2p 2 m 2 /3 Mj 
M ; 16tt 2 / 2 S /* 

The final results for the decay constants are found by combining the wavefunction and vertex cor 
rections: 

Pi+2p 2 \ / 1, 9^ ml/?, , M 2 M? , M 2 

2 m 2 /3 A 2 F(^) A 2 M d . 

35 16. 2 / 2 M 2 -A 2 + 1257 16^ F( X ) I ' (32) 



Pi - 2 P 2 \ ( 1 , # , ml/3 Mj Mj Mj 



VMb^ V M J \ 2 V ^ ; 16^ 2 / 2 b ^ 2 a " 16tt 2 / 2 b M 2 

+ 5 16^ 2 / 2 M 2 -A 2 + 457 16^2 n^)J • (33) 

At this point, let us pause to point out that the nontrivial 1/M dependence enters only only through 
the wavefunction renormalization graphs. The tadpole corrections are only multiplicative — they 
do not see the heavy meson mass. It follows that if the B*Btt couplings are weak, the light and 
heavy quark dependencies of the decay constants will decouple from each other. 

To study the size of 1/M corrections it is useful to consider the ratio U(M) = M/b//b*, which 
is one in the infinite mass limit. In the quenched theory we find 

, N 4« 2 4gg 2 ml/2, Mj 4(972 + 7g 2 ) Mj Mj 

U(M) = 1 + — + -2¥± — Qi — w — <L + yy '- wz> ^— W — i 

K 1 M M 16-7T 2 / 2 p? M 167r 2 / 2 [i 2 

ml/3 A 2 f (^) A 2 M d 

+ 29 16^ 2 / 2 M 2 -A 2 + 891 W^P F{ ~A ); (34) 



G 



while in the full theory 



+ 6 ff 2 (A + 6)*F (£f- s ) 2 ff 2 (A - S) 2 F{^), (35) 

which is easily derived from the work of BG. Here 5 = S s — 5d and in order to keep the expression 
compact we have followed BG in neglecting the pion contribution and approximation to 2 = ^m 2 K , 
but the exact expression is used for numeric work. 



4 Numeric Results 

When studying the error due to quenching, it is interesting to compare our quenched predictions not 
just with the predictions of BG, but also with the predictions of an unquenched two-flavor theory 
with degenerate quark masses. Such a theory should describe most unquenched lattice simulations^ 
and allows one to distinguish between those effects which are due to quenching and those which are 
due to the simplified flavor structure of the lattice simulations. Results for the two flavor theory can 
be extracted from the work of BG by taking into account the SU(3) coefficients. In the quenched 
and two-flavor theories, fs s is obtained from fs by substituting M s = yj2m 2 K — m\ = 680 MeV 
for Md- For convenience, the mass of the light mesons in the quenched and two-flavor theories will 
always be referred to as Md- Since lattice simulations typically have Md > 400 MeV, it will be useful 
to study the Md dependence of our results. 

Before proceeding, it is necessary to specify the various couplings which enter the Lagrangian 
(|l6|). Unfortunately, little information is available. By taking into account chiral loop corrections to 
strong and radiative D* decays, one obtains the constraint 1 15 16) 0.3 < g+{g\—g2)/Mr) < 0.7, with 



the central value being slightly less than 0.5. We will choose g = 0.5 and g\ = g2 as our canonical 
values. QCD sum rules[[l7| |l9|, ^0) and relativistic quark models|2l], ^ favor a smaller value 
of g, g ~ 1/3. They also indicate that g\ — gi is small, on the order of 100 MeV, which is roughly 
consistent with our choice. When necessary, we will assume that both g\ and (72 are of this order. 
This choice is supported by a recent stud y p3| which found 52 — 0.3(? GeV. A lattice study of heavy 
meson decay constants by Bernard et a/.[[24| found p\ + 2p2 = —1.14 GeV, while a similar study by 
the UKQCD collaboration^] found p\ + 2p2 = —0.8 GeV. Because the latter group also studied 
fn* it possible to extract P2 from their data; we find P2 ~ —0.1 GeV, which is consistent with other 
determinations p6| . The remaining couplings 7 and 7^ are unconstrained. Relative to g and gi, they 
are non-leading in 1/N C , so they are likely to be small (in the nucleon system, the singlet coupling 7 
is indeed small[^7|). In view of this, we will take 7 = 7.; — 0, but it should be kept in mind that this 
may be an incorrect assumption. The proper choice of tuq was discussed in Ref. ^| here we choose 
m = 750 MeV. 

We begin by considering the ratio Ri . We first observe that R\ has only one counter-term, whose 
contribution is proportional to (k\i + 2ft3)(m s — md)(l/M[) — I/Mb)- Since we expect from the 
nature of the chiral expansion that the counter-terms will be about the same in the quenched and 
full theories^, it follows that the error in R\ will only weakly depend on the counter-term. With this 
in mind, we neglect it and summarize the results in Table |l||j When the reference values are chosen, 
Ri — 1 is —0.026 in the full theory, —0.053 in the quenched theory, but +0.043 in the two flavor 

3 A more accurate treatment would take into account the difference between the valence and sea quark masses by 
using partially quenched ChPT[ 



4 We do not expect them to - be exactly the same because the two theories have different divergences, but they 
should be close if the lattice is doing a good job of describing QCD 

5 The fact that our result in the three flavor case agrees with that of BG is something of a coincidence, because 
we have chosen to use / = while they use / = fx and we have also retained the pion contribution. 
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Table 1: Quenched and unquenched results for R\. 



Theory 


Ri - 1 




quenched 


-0.2L? 2 


+ 0.20GeV- i g( 5 i-< 72 ) 






+ 0.23<? 7 - 0.097 GeV" 1 ( 5 ( 7 i - 72 )+ 7 ( 5 i - 52 )) 


Nf = 2 


0.17 5 2 


-0.15GeV- 1 g(g 1 -g 2 ) 


N f = 3 


-0.10g 2 


-O.UGeV- 1 g(g 1 -g 2 ) 



Table 2: Numerical results for R\ for various models and couplings. The couplings are chosen as 
follows: Set 1 (BG): g = 0.75, gi = g 2 ; Set 2: g = 0.5, gi = g 2 ; Set 3: g = 0.5, gi - g 2 = 0.2 GeV; 
Set 4: g = 0.5, g x - g 2 = -0.2 GeV. 



Theory 




Hi 


- 1 






Set 1 


Set 2 


Set 3 


Set 4 


quenched 


-0.11 


-0.052 


-0.032 


-0.073 


N f = 2 


+0.087 


+0.043 


+0.028 


+0.058 


N f = 3 


-0.052 


-0.025 


-0.040 


-0.012 



theory. In order to explore the dependence of these results on the couplings, we show in Table g 
the results for several different choices of g and gi. The first set corresponds to the parameters 
chosen by BG, the second set is our canonical one, and the next two explore the dependence on 
gx — g 2 . It can be seen that the error (which is the difference between the full and quenched results) 
is no more than 6% and not half that for the favored choice of couplings. Thus as hoped quenching 
has little effect on R\ — the large tadpole corrections found in Ref. (||) cancel in the ratio and 
the remaining corrections are small. Interestingly, the results in the two-flavor case suggest that 
dynamical simulations will do worse than the quenched simulations. 

The results for R\ can be understood by looking at the ratio fs s //b as a function of M. Figure || 
compares the various theories, neglecting counter-terms. The most striking feature is the "kink" in 



1.4 
1.2 



0.8 
0.6 
0.4 

0.1 0.2 0.3 0.4 0.5 0.6 0.7 O.S 0.9 1 
1/Mb 1/M d 

1/M (GeV 1 ) 

Figure 3: Plot of the ratio fsj Ib in the various theories, 
the three-flavor caseQ, which occurs when A = 5 S , causing the denominator in the argument of one 

6 BG presumably did not comment on this kink because their work only required the ratio at the two physical 
values of M. 
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of the functions F to vanish (c.f. the expression for U(M) in Eq. (|35|)). Physically, this corresponds 
to the case when the B s and the B* are degenerate, leading to an enhancement of f Bs ■ Without this 
kink, the ratio would be a slowly increasing function of M, just as it is in the two-flavor theory. In a 
sense, then, the small error in i?i is a fortunate coicidence. While the details of the curves in Fig. |3| 
depend on our choice of g and gi, the shapes are generic. As g increases, or g\ — gi becomes negative, 
the magnitudes of the slopes increases, so that the curves draw together more at small M without 
changing their shape or relative positions. In particular, there is always a large gap between the 
quenched and full predictions. Thus a general conclusion of this analysis is that quenching decreases 
the ratio /d s //d- In addition, the different M dependence predicted by the three theories makes 
this ratio a test of the chiral description of the decay constants. 

We next consider f B directly. This is primarily for illustrative purposes since the dynami- 
cal theories have an extra counter-term proportional to tr(M + ), making direct comparison to the 
quenched theory difficult; here we will simply neglect the counter-terms. Figure f| then shows both 
Fb = \^Mf B /a and F Bs together with the quenched and two flavor predictions. In the M — ► oo 



1.2 



0.8 
*? 0.6 
0.4 
0.2 


0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 
1/M B 1/M D 

l/M(GeV-') 

Figure 4: Plot of F = vMfs/ot in the various theories, with Md = m^. 

limit the results are consistent with those found in Ref. ||. The behavior of fs is dominated by the 
tree-level term (pi + 2p2)/M. In Fig. ^ we study the Md dependence of the error in f B , comparing 
both the quenched and two flavor theories to the true theory. In the error, the tree-level terms 
cancel, making it easier to see the non-trivial M dependence. One sees that the M and m q depen- 



0.3 

0.2 

ft" 0J 
o 

-0.1 



03 0.4 0.5 0.6 0.7 0.8 0.9 1 
1/M D 

1/M (GeV 1 ) 

Figure 5: Plot of the differences SF B = F B (quenched) - F B (full) and 5F B = F B (N f = 2) — F B (full) 
for various values of Md- 




_i i i i i_ 



1 1 1 1 1 1 


i i i i i i i i i i i i i 

N f =2, M d =m„ 




quenched, M d =m E 




N, = 2, M d =400MeV 




-©-- quenched, M d =400MeV 




-o- N, = 2, M d =600MeV 




-o- quenched, M d =600MeV 




^^^^^^^^^ 








< i . il 



0.1 0.2 
1/M„ 
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dence are closely intertwined, as the both the quenched and two-flavor errors switch sign between 
Md = and Md — 400 MeV. Even at its largest, the quenched error is no more than 7% in the 
region between M — Mb and M — Mo- Also, observe that except for the smallest value of Md, the 
error tends to decrease as M decreases. 

Finally, we examine the 1/M dependence of fg through the quantity U(M). Once again there 
is the problem that the dynamical theories have one more counter-term than the quenched theory. 
This could be eliminated by considering the ratio U s (M)/Ud(M), but we will again simply neglect 
the counter-terms. In Fig. o we show U for various values of the light quark mass. For Md = m^, 





1 1 1 1 1 1 1 1 1 1 1 1 

physical masses 


i i i i i 






Mj = 400. 600 MeV 














N,=3 










quenched 






l<l 


,1 , , , 


N,=2 






0.1 0.2 0.3 0.4 0.5 0.6 0.7 O.S 0.9 
1/M B 1/M D 

1/M (GeV' 1 ) 



Figure 6: Plot of the ratio U(M) — Mfs/fs*- The inset graph shows the M d dependence of the 
quenched and two-flavor theories. As M d increases, the curves of two theories approach each other, 
overlapping near Md = M a and then moving apart again. 

quenching has a pronounced effect. But for larger values of M d , there is almost no difference between 
the quenched and two-flavor theories. Thus lattice simulations are unlikely to observe any effect of 
quenching in U(M). Note that since all the theories must agree at M = oo, the quenching error can 
only increase as M decreases. 

To conclude our discussion of numeric results, let us comment on the uncertainty arising from 
the neglect of the counter-terms. In ChPT for the light mesons, the counter-terms are small, close 
to their natural size of l/(47r) 2 |28| . Goity[Q has observed that choosing /i — 1.5 GeV effectively 
mimics the inclusion of these analytic terms. While we do not know whether this works for the 
heavy mesons, it certainly does generate terms of the natural size. However, it cannot be completely 
correct for the quenched case because of the different behavior of the quenched logarithms under 
scaling. Nonetheless, it serves as a useful way to add counter-terms of the expected size {] Doing 
this, we find that in all the theories R\ is almost independent of /x, changing by only a few percent. 
U{M) also changes very little, again only a few percent, fg and the ratio /bJIb, however, are 
more strongly effected: typical shifts are 10 — 15%. However, the general features of /b s //b that 
were observed in Fig. |§| arc unchanged. 



5 Conclusions 

By generalizing quenched ChPT to include heavy mesons and including 1/M corrections following 
Boyd and GrinsteinjJ], we have studied the 1/M dependence of the quenching errors in lattice 
simulations. We found that unknown counter-terms prevented us from making definitive statements 

7 A more ambitious approach would be to assume that resonances saturate the counter-terms, as they do in the 
light-meson sector |28[ . 
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about the behavior of fs, but it was possible to reduce this uncertainty by taking appropriate ratios, 
such as i?i . Although the long-distance effects of quenching are not independent of the heavy mass 
M, they were fairly insensitive to it. As a result, the error in Ri was found to be small, only a 
few percent. We saw that fsj fs is smaller in the quenched theory, a result that was observed 
earlier in the M — > oo limit ||. We also observed an effect not directly related to quenching, but 
nonetheless missing in lattice simulations, namely the enhancement of fs s when B s and B* are 
nearly degenerate. This enhancement had a strong impact on R\ in the full theory. In future work 
it would be desirable to determine the counter-terms in some way, either from models or directly 
from lattice data. 
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A Renormalized Couplings 



Within the context of dimensional regularization, the singularities of the effective Lagrangian are 
conventionally described in terms of the parameter L(fi) which contains the singularity at D = 4 
(e = 2-D/2, l/e = l/e + log47r+l-7 E ): 



1 



167T 



(A.l) 



With the convention that an arbitrary coupling k is written as k — k r (fi) + kL(/i), the following 
choices render the Lagrangian Eq. ( |l6|) finite. 



kio 
k 3 

~2 



6<?7 , 2 



2^o 
/ 2 



6(571 + 5i7) , 2 



(572 + .927) , 2 



2^o 

P 
2mo 

P 



-2 352 



P 
m|/3 

P 



In addition, it is necessary to add the counter-term 

tr D [H a iv • DbaHb 



3g2 ml/3 



P 

The current Eq. ( |22j ) is renormalized with the choices 

'1 



f>2 



(1 + 3^ )pi +g(3gi -32) 



(1 + 3.9 2 )p 2 +992 



mg/3 



KlO = 697 



Mo 
P' 



(A.2) 
(A.3) 
(A.4) 
(A.5) 
(A.6) 

(A.7) 
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Ml = [5(371 -72) +7(351 -52)] 



K 3 



(.972+752)-^-, 



m|/3 
2^ ' ^ J / 2 



l+o(l + 35 2 )- 



and it is also necessary to rescale a: 

a T (jj,) = a 

Finally, the mass terms in the Lagrangian must also be shifted 

A3 

A 2 = 8.g 



(A.8) 



(A.9) 



I657A2 ■ 2 , 



/ 2 



.2 ™o/3 



/ 2 



(A.10) 
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